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ABSTRACT 

The Large Magellanic Cloud (LMC) star cluster system offers the unique opportunity 
to independently check the accuracy of age (and the corresponding mass) determina- 
tions based on a number of complementary techniques. Using our sophisticated tool 
for star cluster analysis based on broad-band spectral energy distributions (SEDs), 
"AnalySED" , we reanalyse the Hunter et al. (2003) LMC cluster photometry. Our 
main aim is to set the tightest limits yet on the accuracy of absolute age determina- 
tions based on broad-band SEDs, and therefore on the usefulness of such an approach. 
Our broad-band SED fits yield reliable ages, with statistical absolute uncertainties 
within A log(Age/yr) ~ 0.4 overall. The systematic differences we find with respect to 
, previous age determinations are caused by conversions of the observational photometry 

■ to a different filter system, thus leading to systematically inaccurate results. 

"*^h \ The LMC's cluster formation rate (CFR) has been roughly constant outside of 

q h . the well-known age gap between ~ 3 and 13 Gyr, when the CFR was a factor of 

1 1 ^ 5 lower. Using a simple approach to derive the characteristic cluster disruption 

time-scale, we find that log(ifVyr) = 9.9 ± 0.1, where t dis = tf 8 (M cl /lO 4 M ) - 62 . 
This long characteristic disruption time-scale implies that we are observing the initial 
cluster mass function (CMF). We conclude that while the older cluster (sub)samples 
show CMF slopes that are fully consistent with the a ~ —2 slopes generally observed 
in young star cluster systems, the youngest mass and luminosity-limited LMC cluster 
subsets show shallower slopes (at least below masses of a few xlO 3 M ), which is 
contrary to dynamical expectations. This may imply that the initial CMF slope of the 
LMC cluster system as a whole is not well represented by a power-law, although we 
cannot disentangle the unbound from the bound clusters at the youngest ages. 

Key words: methods: data analysis, Magellanic Clouds, galaxies: star clusters, galax- 
ies: stellar content 
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1 INTRODUCTION 

The debate regarding the true underlying mass function of 
young star cluster systems is presently very much alive, be- 
cause it bears on the very essence of the star-forming pro- 
cess, as well as on the formation, assembly history and evo- 
lution of the clusters' host galaxies over cosmic time. Al- 
though observational evidence of young cluster systems in 
many star cluster-forming interacting and starburst galaxies 
appears to indicate that they are well represented by power- 
law cluster luminosity functions (CLFs; see de Grijs et al. 
2003d for a recent comprehensive review), rapid changes to 
the stellar population properties, combined with a (possi- 
bly significant) age range within a given cluster system may 
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conspire to give rise to a power-law-like CLF, even if the 
true underlying cluster mass function (CMF) is not a power 
law (Meurer 1995; Miller et al. 1997; Fritze-v. Alvensleben 
1998, 1999; de Grijs et al. 2001, 2003a,b; Hunter et al. 2003, 
hereafter H03). It is, therefore, obviously very important to 
age date the individual clusters and to correct the observed 
CLF to a common age, before interpreting the cluster lu- 
minosities in terms of the corresponding mass distribution 
(Fritze-v. Alvensleben 1999; de Grijs et al. 2001, 2003a,b; 
H03). This is particularly important if the age spread within 
a cluster system is a significant fraction of the system's mean 
age. 

Even with the high spatial resolution offered by the ad- 
vanced suite of instruments onboard the Hubble Space Tele- 
scope, extragalactic star clusters at distances beyond the 
Local Group can only be studied through their integrated 
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properties. However, we do have access to a statistically 
significant young to intermediate-age resolved star cluster 
system on our doorstep, notably in the Large Magellanic 
Cloud (LMC). By combining integrated properties with re- 
solved stellar population studies, the LMC cluster system 
offers the unique chance to independently check the accu- 
racy of age (and corresponding mass) determinations based 
on broad-band spectral energy distributions (SEDs) . This is 
what we set out to do in this paper, based on the largest, 
most homogeneous and complete set of integrated LMC clus- 
ter photometry presently available. With our new cluster age 
estimates, spanning ages from a few Myr up to 10 Gyr, de- 
termined in an internally consistent fashion and with a firm 
handle on the associated uncertainties, we now also have the 
largest, homogeneous sample of cluster masses with well- 
defined uncertainties to date, which we use to explore the 
evolution of the CMF. 

In Section 2 we discuss the details of the LMC cluster 
sample we base our analysis on, and explore the uncertain- 
ties associated with the age and mass determinations. Sec- 
tion 3 discusses the formation and disruption history of the 
LMC star cluster sample. We use this as our basis for the 
interpretation of the LMC cluster mass distribution in terms 
of the initial and evolved distributions of cluster masses in 
Section 4. Finally, in Section 5 we summarise our results and 
conclusions. 



2 DATA AND ACCURACY ANALYSIS 

The basis for our detailed comparison of the age estimates 
for the LMC star cluster system is provided by (i) the 
UBVR broad-band SEDs of H03, based on Massey's (2002) 
CCD survey of the Magellanic Clouds, and (ii) the OGLE-II 
data set (Udalski, Kubiak & Szymahski 1997), and in partic- 
ular their homogeneously determined ages for some 600 clus- 
ters, using colour-magnitude diagrams (CMDs; Pietrzyhski 
& Udalski 2000). 

H03 determined the ages for the individual LMC clus- 
ters by comparing their broad-band SEDs with a variety of 
cluster evolutionary models. Specifically, they used the Star- 
burst99 (Leitherer et al. 1999) with Z = 0.008 for ages up 
to 1 Gyr, while for older ages they used a combination of 
the UBV colours of Searle, Sargent & Bagnuolo (1973), the 
(V — R)c colours for globular clusters from Reed (1985), 
and the Chariot & Bruzual (1991) simple stellar population 
(SSP) models for the evolution of the U-band luminosity in 
the age range from 1 to 10 Gyr. They then used a combi- 
nation of two colour-colour diagrams, (U — B) vs. (B — V) 
and (B — V) vs. (V — R), to obtain their age estimates, by 
comparing the cluster positions to the model tracks in those 
diagrams. 

The colours and magnitudes of their sample clusters 
were corrected for reddening using a blanket reddening of 
E(B—V) = 0.13 mag for all clusters, where they adopted the 
extinction curve of Cardelli et al. [1989; Rv = Ay/E(B — 
10=3.10]. 

The Pietrzyhski & Udalski (2000) ages are based on 
detailed isochrone fits to the OGLE-II colour-magnitude 
data, using the Padova isochrones for Z — 0.008. They 
dereddened their photometry using the extinction values 
for 84 LMC subfields determined by Udalski et al. (1999), 



E(B - V) = 0.143 ± 0.020 mag. These extinction values 
are based on Schlegel et al.'s (1998) reddening maps for the 
Galactic foreground extinction, on average E(B — V) = 0.075 
mag, and their Rv = 3.24. 

Finally, here we reanalyse the Massey (2002) data set, 
also assuming Z = 0.008, for comparison purposes (although 
our analysis tool allows us to leave the metallicity as a free 
parameter; see Section 2.2). For the total extinction towards 
the LMC clusters, we assumed E(B — V) = 0.10 mag, using 
the Calzetti attenuation law (Calzetti 1997, 2001; Calzetti 
et al. 2000; Leitherer et al. 2002) with R v = 4.05. This 
corresponds to E(B — V) ~ 0.13 mag for both the Cardelli 
et al. (1989) and the Schlegel et al. (1998) extinction laws. 

In a series of recent papers, we developed a sophisticated 
tool for star cluster analysis based on broad-band SEDs, 
"AnalySED", which we tested extensively both internally 
(de Grijs et al. 2003c,d; Anders et al. 2004) and externally 
(de Grijs et al. 2005), using both theoretical and observed 
young to intermediate-age (t < 3 x 10 yr) star cluster SEDs, 
and the GALEV SSP models (Kurth et al. 1999; Schulz et al. 
2002). We increased the accuracy for younger ages by the 
inclusion of an extensive set of nebular emission lines, as 
well as gaseous continuum emission (Anders & Fritze-v. Al- 
vensleben 2003). We concluded that the relative ages within 
a given cluster system can be determined to a very high ac- 
curacy depending on the specific combination of passbands 
used (Anders et al. 2004) . Even when comparing the results 
of different groups using the same data set, we can retrieve 
prominent features in the cluster age distribution to within 
a t = A(log(Age/yr)) < 0.35 (de Grijs et al. 2005), which 
confirms that we understand the uncertainties associated 
with the use of our AnalySED tool to a very high degree. 

The aim of this section is to compare our new age deter- 
minations with (i) those of H03 using the same data set but 
a different approach, and (ii) those of Pietrzyhski & Udalski 
(2000) using the independent (and presumably more accu- 
rate) method of CMD fitting. This will allow us to set the 
tightest limits yet on the accuracy of absolute age determi- 
nations based on broad-band SEDs, and therefore on the 
usefulness of such an approach in general. 

First, we compare the independently determined ages of 
the clusters in common between the H03 and OGLE-II re- 
sults. Fig. la shows the distribution of the data points in the 
relevant parameter space, where the dashed line represents 
equality between both age estimates for a given cluster. It is 
immediately clear that there is a systematic difference be- 
tween the H03 and OGLE-II ages, the magnitude of which 
is a clear function of age. The effect is indeed significant; 
a straightforward linear regression results in a slope that is 
significantly different from unity, 1.47±0.07. We will explore 
the origin of this systematic effect in Section 2.1 below. It is 
also clear that, in units of log(Age/yr), the intrinsic scatter 
in the data points is somewhat greater at younger ages; this 
is merely an effect of the poorer intrinsic age accuracy for 
older ages imposed by the models. 

Secondly, we compare our redetermined ages to both 
those of OGLE-II and H03, in Figs, lb and c, respectively. 
We point out that because of the limiting U-band magnitude 
of the OGLE-II data, their age estimates are only reliable 
for ages < 1 Gyr. Consequently, the comparison between the 
OGLE-II CMD-based ages and our age determinations in 
this paper are only valid for ages < 1 Gyr, as is also evident 
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Figure 1. Comparison of age determination for LMC clusters in common between Hunter ct al. (2003), the OGLE-II team of Pietrzynski 
& Udalski (2000), and this paper. The assumptions on the overall metallicity and extinction properties of the cluster sample arc indicated 
in the legend. 



from the lack of data points for log(Age/yr) > 9.0 in Fig. 
lb. For the older sample clusters, for which age estimates 
are available based on spectral features or CMD analysis, 
the comparison is extended on a one-by-one basis in Section 
2.3. 

The systematic effect shown by the H03 data seems to 
have disappered when we compare our results to those of the 
OGLE-II team; the resulting slope is 1.05 ±0.05, with an x- 
axis intercept at log(Age/yr) = —0.11 ± 0.39. This excellent 
agreement indicates that our assumption that the overall ex- 
tinction is well represented by E(B — V) = 0.10 mag (assum- 
ing the Calzetti attenuation law) is well justified; a change 
in extinction of \AE(B — V)\ = 0.05 would correspond to 
an offset from the line of equality of |A log(Age/yr)| ~ 0.10. 
This effect is a direct result of the age-extinction, and possi- 
bly (to a lesser extent) also of the age-metallicity degeneracy 
(see Section 2.2). The systematic effect with respect to the 
H03 results is still visible, although the scatter in the results 



is - perhaps not surprisingly because we used the same data 
set - much smaller than when we compare our results to 
those of OGLE-II. Table 1 contains the full numerical details 
for a quantitative comparison among the age determinations 
based on the samples discussed here. We note that any re- 
maining offset between our results and those of the OGLE-II 
team, despite having used the same assumptions for the total 
extinction towards the clusters, may have been introduced 
by the exact value for the distance modulus assumed by 
Pietrzynski & Udalski (2000), m- M = 18.23 mag, which 
is an essential ingredient for age determinations based on 
CMD fitting. If they had assumed a nominal m — M = 18.50 
mag, their resulting age estimates would have been younger 
by Alog(Age/yr) ~ 0.2 — 0.4, depending on the age range 
considered (larger offsets result for younger ages). 
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Table 1. Quantitative comparison among age determinations. 
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Figure 2. (a)-(c) The effects of colour transformations between filter systems on the derived ages for the Massey (2002) LMC cluster 
photometry, assuming Z = 0.008 and E(B — V) = 0.10 mag. We have chosen to omit the formal uncertainties on the data points for 
reasons of clarity, (d) Comparison of the corrected age determinations for the H03 sample to our redeterminations in this paper. 



2.1 Filter transmission curves 

In the previous section we noted a significant systematic ef- 
fect between the age differences of H03 on the one hand, and 
those of both the OGLE-II team and our own redetermina- 
tions on the other. H03 converted the Starburst99 Johnson 
(V — R) colour to the Cousins system. As we showed in de 
Grijs et al. (2005; their figs. 10 and 11), filter conversions 



may be responsible for significant differences in the result- 
ing age determinations. In order to investigate the possibil- 
ity of this effect playing a role for these data, we carefully 
analysed the properties of the original and converted sys- 
tems used. Massey (2002) calibrated his photometry using 
the Landolt (1992) standard stars, so that his photometry 
should be analysed using the filter set used to obtain the 
standard-star photometry. The Landolt (1992) filter curves 
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are very close to the filter transmission curves used by Holtz- 
man et al. (1995), who also kindly supplied us with the orig- 
inal, unpublished Landolt KPNO curves (J. Holtzman, priv. 
comm.). For completeness, we will therefore assess the dif- 
ferences in the resulting ages for our LMC cluster sample 
based on using the "standard" Johnson-Cousins system (as 
done by H03), the curves used by Holtzman et al. (1995), 
and the original KPNO curves used by Landolt (1992). 

The results of this are shown in Fig. 2. Fig. 2c shows 
that the filter transmission curves used by Holtzman et al. 
(1995) and Landolt (1992), respectively, indeed yield rela- 
tively similar age estimates. The most important compar- 
ison figure is displayed in panel b, however. The system- 
atic trend seen here mimics that seen in Figs, la and c, in 
the sense that if one uses the Johnson- Cousins filter system 
(even if based on the appropriate conversion equations) in- 
stead of the native Landolt KPNO system, one will obtain 
lower ages than expected at the young-age end of the age 
range covered by our LMC cluster sample, and vice versa. 
To provide further evidence for this scenario, we applied a 
simple linear regression to the data points in Fig. 2c (shown 
by the dotted line), and applied the resulting equation as 
a first-order correction to H03's age determinations in both 
Figs, la and c. The resulting, corrected age comparison is 
shown in Fig. 2d, while the numerical values are once again 
included in Table 1. It thus appears that the systematic dif- 
ferences in LMC cluster ages between the H03 results and 
those of OGLE-II and ourselves are indeed caused by their 
conversions of the photometry to a different filter system. 

2.2 The age-metallicity and age-extinction 
degeneracies 

In the discussion of our results, we have thus far assumed 
a fixed metallicity of Z — 0.008 and a fixed extinction of 
E(B — V) = 0.10 mag (assuming the Calzetti attenuation 
law) for all of our sample clusters. The AnalySED tool has 
been developed to also provide independent information on 
a cluster's metallicity and extinction, provided that a signif- 
icant number of data points defining the SED are available 
to match the number of free parameters (cf. de Grijs et al. 
2003c,d; Anders et al. 2004). 

Therefore, we have also reanalysed the LMC cluster 
photometry assuming (i) variable metallicity and E(B — 
V) = 0.10 mag, (ii) Z = 0.008 and variable extinction, 
and (iii) variable extinction and metallicity, the results of 
which are shown in Figs. 3a, b and c, respectively. From close 
scrutiny of these panels, it appears that the age-metallicity 
degeneracy is least important for this cluster sample; this 
implies that the assumption of Z = 0.008 adopted by H03, 
OGLE-II and ourselves in the previous sections is a reason- 
able approximation of the average LMC cluster metallicity. 
We also note that adopting Z — 0.008 as a reasonable ap- 
proximation for the average cluster metallicity is supported 
- at least for clusters younger than ~ 3 Gyr - by a number of 
spectroscopic studies based on individual cluster stars (see, 
e.g., Olszewski et al. 1991, their fig. 11). 

The effects of the age-extinction degeneracy are dra- 
matic for this sample, however. On average, we find that if 
we had not assumed a fixed extinction value, the resulting 
ages would have been smaller (although there remains sig- 
nificant scatter about the dashed line of equality). However, 



we caution that some of this scatter is likely artificial; as we 
showed in de Grijs et al. (2003d), the age-extinction degen- 
eracy is artificially worsened (in the sense that clusters are 
artificially assigned younger ages) if the wavelength base- 
line used for the multi-parameter cluster SED analysis lacks 
wavelength coverage longward of the R filter, as is the case 
here. In essence, this is shown by the enhanced population 
of the two youngest age bins (at 4 and 8 Myr) in Figs. 3b 
and c for the models without restrictions on E(B — V). This 
situation can be avoided successfully by adopting a fixed 
value for either the cluster population's (average) extinction 
or metallicity (de Grijs et al. 2003d). 

Thus, in view of these considerations, this implies that 
our results are indeed very comparable to those obtained 
from the OGLE-II data. 

2.3 Other comparisons 

Finally, to assess the robustness of the individual ages, we 
compare the LMC cluster age determinations of H03 (un- 
corrected), OGLE-II and ourselves with those obtained for 
smaller cluster samples using a variety of techniques. In Fig. 
4 we have included those samples containing > 5 clusters in 
common with our sample for which the authors determined 
their ages based on broad-band photometry. With the ex- 
ception of the Meurer et al. (1990) results, who used the 
ultraviolet colours C(18 — 31)o and C(25 — 31)o, all other 
age determinations are based on blue optical colours. Bica 
et al. (1986, 1990) used the C(43 - 45) vs. C(45 - H/3) 
vs. SWB type (Searle, Wilkinson & Bagnuolo 1980) diag- 
nostics, while Chiosi et al. (1988; UBV clours) and Santos 
et al. (1995; (17 — B)o colours) limited themselves to the 
"standard" set of broad-band filters. 

Although the Meurer et al. (1990) ages are clearly sys- 
tematically greater than any of the ages determined by H03, 
OGLE-II or us, the other age determinations are statisti- 
cally similar to those discussed elsewhere in this paper. This 
is encouraging regarding the robustness of our results; the 
relatively large scatter is a side effect of the poor age reso- 
lution for broad-band colours only. The fact that Meurer et 
al.'s (1990) ultraviolet-based age determinations are clearly 
greater than ours underscores the difficulties involved in 
using this wavelength range alone for age determinations. 
The main uncertainties inherent to this at ultraviolet wave- 
lengths are the uncertain stellar population synthesis and 
uncertain extinction corrections. 

In Fig. 5 we show a similar comparison between the ages 
discussed elsewhere in this paper, and those obtained inde- 
pendently based on spectral features or CMD fits. The simi- 
larities between our and H03's broad-band SED fits and the 
OGLE-II fits on the one hand, and the determinations from 
the literature are striking, as is the relatively small scatter 
about the (dashed) lines of equality. Both the simplest of 
these age determinations, Elson & Fall's (1988) ages based 
on the clusters' main-sequence turn-off magnitudes (assum- 
ing a common distance to their entire sample), as well as 
the CMD fits by Olszewski et al. (1991) and Vallenari et 
al. (1998), and the more sophisticated determinations based 
on distinct spectral features by Beasley et al. (2002; see fig- 
ure legend) and Santos & Piatti (2004; based on equivalent 
width measurements of the integrated spectra) match our 
broad-band SED ages very well. We also emphasize that 
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Figure 3. The effects of (a) the age-metallicity, (b) the age-extinction, and (c) the combined degeneracies. Again, error bars have been 
omitted for reasons of clarity. 



Barbaro & Olivi's (1991) ages based on ultraviolet spectra 
return reliable ages, as opposed to those based on the ultra- 
violet colours discussed in relation to Fig. 4. H03, based on 
a comparison involving fewer objects taken from the litera- 
ture, also concluded that the general match between their 
results and CMD-based ages was satisfactory, although they 
stated that for the oldest cluster colours tend to underesti- 
mate the ages compared to CMD fits. However, we have now 
shown that most (but possibly not all; see below) of this dis- 
crepancy was most likely caused by the systematic effects in- 
troduced by the filter conversion they employed. The small 
number of clusters available for this comparison (and the 
smaller number used by H03) did not allow H03 to notice the 
systematic offset discussed in the previous sections. Based 
on the small number of clusters aged ~ 1 Gyr for which 
comparison data is available in the literature, we tentatively 
conclude that our broad-band SED analysis may have led 
to underestimates of the cluster ages for log(Age/yr) ~ 9, 



in a similar sense as seen by H03. This tentative conclusion 
is supported by Fig. 5a, where our new age determinations 
appear somewhat lower than those determined from CMD 
analysis or spectral features for these ages. 



Based on this comparison with independently deter- 
mined ages for subsamples of LMC clusters, we conclude 
that our broad-band SED fits yield reliable ages, with statis- 
tical absolute uncertainties within A log(Agc/yr) ~ 0.4 over- 
all, based on the scatter about the lines of equality. Thus, 
in addition to our conclusion in de Grijs et al. (2005) that 
we can retrieve prominent features in the cluster age distri- 
bution to within A(log(Age/yr)) < 0.35, we have now also 
shown that the intrinsic statistical uncertainties involved in 
cluster age determinations based on broad-band SEDs are 
of a very similar magnitude. 
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Figure 4. Comparison with other age determinations, based on broad-band photometry. 



3 CLUSTER FORMATION AND DISRUPTION 
RATES 

Using the newly determined and improved age estimates for 
the largest LMC cluster sample spanning the most extensive 
age and mass ranges to date, we now have the means to 
constrain the past cluster formation rate (CFR), as well as 
the effects of cluster disruption, to unprecedented detail. 

3.1 Derivation of the characteristic cluster 
disruption time-scale 

In Fig. 6a we display the number of clusters formed per 
unit age range. Evolutionary fading, as predicted by stel- 
lar population synthesis models, will cause this number to 
slowly decline from the youngest ages upward. This effect is 
shown by the dashed line in Fig. 6a, of which the slope is 
entirely determined by the details of stellar population syn- 
thesis; we have only shifted this line to best match the data 
points for log(Age/yr) < 8. Both internal and external ef- 



fects, such as two-body relaxation, disk and bulge shocking, 
and the tidal effects caused by the underlying galactic grav- 
itational potential (even in the low-density environment of 
the LMC; cf. Lamers, Gieles & Portegies Zwart 2005), lead- 
ing to tidal stripping and to evaporation of a fraction of the 
low-mass cluster stars, will result in the (gradual) dissolu- 
tion of star clusters. Simple (instantaneous) disruption the- 
ory (Boutloukos & Lamers 2003) predicts that, assuming a 
constant CFR and that the characteristic cluster disruption 
time-scale is a function of cluster mass, the effects of cluster 
disruption will dominate from a certain age onwards, giving 
rise to the characteristic double power-law seen in Fig. 6a. 

The location of the crossing point of the fading and 
disruption lines in Fig. 6a, denoted by "tcross", is a key in- 
gredient to derive the characteristic cluster disruption time- 
scale for a given cluster population (Boutloukos & Lamers 
2003). The other parameters required to derive this time- 
scale are the limiting magnitude, My = —3.5 mag (H03), 
E(B - V) = 0.10 mag or A v = 0.41 mag (Section 2), 
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Figure 5. Comparison with other age determinations, based on either spectral features or CMD fitting. 



(rn — M)lmo = 18.48 mag, and the reference magnitude, 
My ! = -11.3668 mag, for a 10 4 M Q cluster at an age 
of 10 8 yr, which we obtained from the GALEV models for 
Z — 0.4 Z , assuming a "standard" Salpeter-like stellar 
IMF from 0.15 to ~ 70 M©. For the crossing time we find 
log(t cross /yr) = 8.10 ±0.05, where the uncertainty is a com- 
bination of the uncertainty in the exact vertical level of the 
fading line, and the uncertainty regarding exactly when clus- 
ter disruption is significantly more important than fading. 
In order to assess the latter uncertainty, we applied a linear 
regression to the data points for log(Age/yr) > 8.2 and 8.4, 
respectively, resulting in the dotted and dash-dotted fits in 
Fig. 6a. We note that for these fits we have omitted the data 
points in the LMC's age gap, between 3 and 13 Gyr (see Sec- 
tion 3.2), since these would obviously violate the underlying 
assumption of a constant CFR. 

In addition, we need the slope of the fading line (—0.648; 
from stellar population synthesis modelling) and the exact 
mass dependence of the disruption time-scale, 7, defined 



as tdis oc MJ, (Boutloukos & Lamers 2003), in order to 
derive the characteristic disruption time-scale. Boutloukos 
& Lamers (2003) derived 7 = 0.62 from observations of a 
small but diverse sample of galaxies containing rich clus- 
ter systems. Using these parameters as our input, we de- 
rive a characteristic cluster disruption time-scale for a 10 4 
M cluster of log(tf B /yr) = 9.9 ± 0.1. This is in good 
agreement with Boutloukos & Lamers (2002), who found 
log(£4 i 7y r ) = 9- 7 ± 0.3 for a smaller sample of 478 clusters 
within 5 kpc from the centre of the LMC, in the age range 
7.8 < log(Agc/yr) < 10.0 (data from Hodge 1988), while our 
result is also qualitatively consistent with H03 who noticed 
very little destruction of clusters on the high-mass end. This 
is not unexpected, considering the low-density environment 
in which these clusters are found (cf. Lamers et al. 2005). 

Even with recent improvements (e.g., Gieles et al. 2005) 
to the simple model of Boutloukos & Lamers (2003), the 
characteristic cluster destruction time-scales resulting from 
a more sophisticated, non-instantaneous destruction process 
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oldest clusters that may possibly have formed at a different 
CFR (if we leave out the oldest clusters, the main effect is 
that the slope of the disruption line becomes less well con- 
strained); (ii) the slope of the disruption line in either case 
corresponds to a mass dependence of the LMC cluster dis- 



4 5 
log( M a / M e ) 



Figure 6. (a) - The LMC cluster formation rate (in number of 
clusters per Myr) as a function of age. The dashed line is the 
least-squares power-law fit to the fading, non-disrupted clusters, 
for a constant ongoing cluster formation rate. The dotted and 
dash-dotted lines are the disruption lines for the most likely age 
ranges where disruption may dominate evolutionary fading, as 
discussed in the text; (b) — Mass spectrum of the LMC clusters 
(number of clusters per unit mass range); the dotted line rep- 
resents a power-law fit to the CMF that is as yet unaffected by 
disruption (see text); the slope of the dash-dotted line (which was 
shifted vertically to match the observational data) is entirely de- 
termined by the initial CMF slope determined from the dotted 
line, and stellar population synthesis, as described in the text. 



are very similar to those based on the simple method (see 
Lamers et al. 2005 for a detailed comparison, in particular 
their table 1). 

Finally, for completeness and a full assessment of the 
uncertainties involved, we need to address the effects of 
possibly having underestimated the ages of the oldest sam- 
ple clusters. The ages of the oldest LMC clusters are well- 
known from detailed Hubble Space Telescope-based stud- 
ies (e.g., Olsen et al. 1998) to all have essentially the same 
age of log(Age/yr) ~ 10.1. Therefore we reanalyzed Fig. 
6b, now assuming this uniformly old age for all clusters 
for which we determined log(Age/yr) > 9.5. Under this as- 
sumption, the resulting crossing time requires adjustment by 
A log (t cr oss /yr) — —0.05, leading to a lower limit to the char- 
acteristic disruption time-scale of log(£4 ls /yr) = 9.8 ± 0.1, 
which is still in very good agreement with the earlier deter- 
mination by Boutloukos & Lamers (2002). 

Our analysis in this section depends on the key assump- 
tion that the CFR has remained approximately constant for 
the entire lifetime of the LMC star cluster system, with the 
exception of the time span covered by the age gap. That 
this assumption is justified to first order is supported by 
two additional lines of evidence: (i) the disruption lines in 
Fig. 6a are representative of the full LMC cluster sample for 
ages of log(Age/yr) > 8.1, whether or not we include the 
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7=0.56±0.07 



which is - within 



ruption time-scale of ta ls ^ j« c1 
the uncertainties - well inside the approximately universal 
proportionality found by Boutloukos & Lamers (2003) for a 
range of galaxies in the local Universe, 7 = 0.62 ± 0.06. 



3.2 The LMC cluster formation rate 

In Fig. 7 we show the distribution of the full LMC cluster 
sample in the (age vs. mass) plane. For reasons of clarity, 
we have not included error bars in this figure. Overplot- 
ted is the 50 per cent completeness limit based on stellar 
population synthesis for single- burst ("simple") stellar pop- 
ulations (appropriate in the context of star clusters), for an 
estimated ~ 50 per cent completeness limit of My = —3.5 
mag (H03), assuming a distance modulus to the LMC of 
(m — M)o = 18.48 mag, and no extinction. For a nominal 
extinction of Ay =0.1 mag (assuming the Calzetti extinc- 
tion law), the equivalent detection limit will shift to higher 
masses by Alog(M c i/M ) = 0.04, which is well within the 
uncertainties associated with our mass determinations. 

For a constant CFR, we expect the number of data 
points to increase gradually and evenly from young to old 
ages in the logarithmic representation of Fig. 7. Any obvi- 
ously clumped subsets of data points may indicate either 
artifacts caused by our computational routines, or real devi- 
ations from a constant CFR. The effects of using interpola- 
tion of discrete isochrones for the determinations of cluster 
ages, and the resulting artifacts in age-mass space have been 
discussed extensively by Bastian et al. (2005) and Gieles 
et al. (2005). Based on their discussions, two features in 
Fig. 7 in particular can be attributed to artifacts caused 
by the fitting procedures: (i) the large density of clusters 
at log(Age/yr) = 6.6 is simply caused by the fact that our 
youngest isochrone is at an age of 4 Myr (we are limited by 
the age range spanned by the Padova isochrones on which 
the GALEV SSP models are based), and younger clusters 
would therefore be assigned the minimum model age; (ii) the 
apparent overdensity around log(Age/yr) = 7.2 is caused by 
the discreteness of the isochrones in this age range, where 
rapid changes occur in realistic stellar populations; it was 
also noted by Gieles et al. (2005) and Fall et al. (2005). For- 
tunately, however, from their age and mass determinations 
of artificial cluster populations using very similar procedures 
as done here, Gieles et al. (2005) find that they can retrieve 
the ages and masses of, respectively, 87 and 97 per cent of 
their input clusters to within 0.4 dex in log(Age/yr) and 
log(Mci/M Q ), respectively (see also Fall et al. 2005). These 
retrieval rates are based on the assumption of a constant 
CFR, and an initial power-law CMF with exponent a = —2. 

In addition to these artifacts, we inspected two other 
apparent overdensities in Fig. 7, at (7.8 < log(Age/yr) < 
8.0,2.8 < log(M c i/M ) < 3.4) and 8.2 < log(Age/yr) < 8.4 
(for the full mass range) , by displaying the locations of these 
cluster subsamples across the LMC. Neither of these ap- 
parent overdensities are associated with localised enhanced 
CFRs. In fact, any sufficiently large age and mass range 
(i.e., covering at least twice the typical uncertainty in age 
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and mass, following the Nyquist sampling theorem) shows 
a fairly homogeneous and relatively continuous (in terms of 
the cluster ages) distribution of star clusters across the disk 
of the LMC. We note, as a caveat, that the H03 cluster sam- 
ple does not cover the entire LMC (e.g., it does not cover 
the dense LMC bar region), although Massey (2002) and 
H03 argue that their fields are representative of the LMC 
stellar (and cluster) population as a whole. 

The last tidal encounter between the LMC and the 
Small Magellanic Cloud (SMC) occurred most likely about 
0.2 to 0.5 Gyr ago (e.g., Heller & Rohlfs 1994; Gardiner 
& Noguchi 1996). One rotation period of the LMC cor- 
responds to ~ 250 Myr (Grebel & Brandner 1999), or 
log(tLMc/y r ) — 8-4- This implies that if an apparent fea- 
ture in age-mass space at an age of log(Age/yr) ~ 7.9 was 
triggered by the LMC-SMC encounter, but corresponds to a 
spatial distribution of its clusters that is smooth and homo- 
geneous across the LMC disk (i.e., at all radii), rotational 
mixing will not have had enough time for this to have hap- 
pened, or at best marginally so. While such a gravitational 
interaction could, in principle, have triggered a galaxy-wide 
burst of cluster formation in the LMC, we would still ex- 
pect to see a clearly enhanced population of clusters in a 
relatively narrow age range in such a case (cf. de Grijs et al. 
2003a,b for the equivalent diagnostics in M82). Lacking bet- 
ter age resolution for a larger distinct cluster subpopulation, 
we cannot conclude with sufficient certainty that the last 
encounter between the LMC and the SMC was significantly 
strong to have triggered the formation of a well-defined clus- 
ter population in the LMC as a result; alternative triggering 
mechanisms may have been responsible for the formation of 
the smoothly distributed clusters at young(er) ages. 

The only real unrferdensity of data points is the well- 
known LMC cluster age gap, between ~ 3 and 13 Gyr (e.g., 
Da Costa 1991; Geisler et al. 1997; Rich et al. 2001; Piatti 
et al. 2002; Bekki et al. 2004), which is well reproduced in 
our broad-band age redeterminations. We note that possible 
other underdensities in Fig. 7 occur (i) close to the 50 per 
cent completeness limit and (ii) for low masses, where we 
have to take stochastic sampling effects of the stellar initial 
mass function (IMF) into account; these effects act in the 
sense that a cluster of low mass may be dominated by a small 
number of high(er)-mass, and therefore high(er)-luminosity 
stars, which in turn will cause the photometric age to be 
(somewhat) overestimated. The exact amount of this de- 
pends sensitively on the cluster mass (see Section 4 below) 
and age. In particular, stochastic effects will be most pro- 
nounced in the age ranges where rapid evolutionary stages 
show up in stellar models, such as the appearance of red 
supergiants at ages around 8 — 12 Myr (e.g., Girardi & Bica 
1993; Santos & Frogel 1997). Clusters that have reached ages 
of several x 10 7 yr will be populated by a significant fraction 
of post-main sequence stars, and suffer less from stochastic 
effects in optical passbands, even if they are of low mass 
(e.g., Girardi & Bica 1993; Santos & Frogel 1997). The ap- 
pearance of thermally-pulsing asymptotic giant branch stars 
around ages of ~ 1 Gyr can potentially give rise to large 
stochastic IMF sampling effects. However, these are most 
important at near-infrared wavelengths, and of lower im- 
portance at the optical wavelengths used for the age and 
mass estimates in this paper. Additional, yet circumstan- 
tial supporting evidence for this statement is provided by 
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Figure 7. Distribution of the LMC clusters in the (age vs. mass) 
plane. Overplotted is the expected detection limit based on stel- 
lar population synthesis for a completeness limit of My = —3.5 
mag (H03), assuming no extinction. For a nominal extinction of 
Ay = 0.1 mag (assuming the Calzetti attenuation law), the detec- 
tion limit will shift to higher masses by Alog(M c i/M0) = 0.04, 
which is well within the uncertainties associated with our mass 
determinations. The features around 10 Myr are caused by the 
appearance of red supergiants in the models. 



the age and mass distribution shown in Fig. 7, where we do 
not detect any significant under or overdensity in mass of 
clusters at or close to ages of log(Age/yr) ~ 9.0, although 
such effects could - to some extent - be masked by the rel- 
atively small number (< 100) of clusters in the logarithmic 
age interval [8.8,9.2]. 

An alternative method to portray the CFR is shown 
in Fig. 6a, where we display the number of clusters formed 
per unit age range. If we omit the three data points in the 
age gap, the remaining data points for log(Agr/yr) > 8 are 
well approximated by a linear relationship. Since a linear 
relationship is predicted if the CFR has been constant, this 
approach shows the approximate validity of this assumption 
once again. If we assume that the CFR in the cluster age 
gap has also been constant, although at a (much) lower level 
(cf. H03), then we derive that the ratio of the CFR before 
and after the cluster age gap to that in the period between 
~ 3 and 13 Gyr, is about (5 ± 1) : 1. H03, considering the 
upper cluster masses, and assuming that the mass of the 
most massive cluster is determined by size-of-sample effects 
(which is approximately valid; H03), derived an equivalent 
ratio of ~ 10, although with large uncertainties. In view 
of the independent methods used by H03 and us, based on 
entirely different diagnostics, these results are in reasonable 
agreement. 

Finally, H03 suggested that the oldest sample clusters 
represent a separate population from the younger LMC clus- 
ters. They state that lower-mass old clusters could have been 
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observed but are not in practice. However, our new age and 
mass determinations do not concur with this result; Fig. 6 
shows that our 50 per cent completeness limit describes the 
lower mass limit of the entire cluster sample up to the old- 
est ages very well. Similarly, the upper mass limit of the 
oldest clusters is commensurate with the upper mass limit 
expected from size-of-sample effects (cf. H03), and does not 
require that the oldest clusters represent a separate popula- 
tion. However, because of the small number of old clusters, 
we cannot draw any firmer conclusions on this issue based 
on the information at hand. 



4 IMPLICATIONS FOR THE CLUSTER MASS 
FUNCTION AS A FUNCTION OF AGE 

In Fig. 6b we show the LMC cluster mass distribution in 
a way that allows us to assess the importance of disruption 
processes for the sample, in the presence of an age-dependent 
detection (completeness) limit (see Boutloukos & Lamers 
2003). As we showed in de Grijs et al. (2003d), the slope 
of the distribution for the highest masses is in essence a 
projection of the initial CMF, provided that we can prove 
that (semi-instantaneous) disruption will not yet have had 
the time to act on these masses. 

In Table 2 we list the derived slopes for mass ranges 
from the minimum masses indicated to the highest masses 
present in our sample. We note that, for log(M c i/Mo) > 3.0 
the CMF slopes derived are very stable (except for the 
highest mass ranges, where the small number of clusters 
affects the results), resulting in an intial CMF slope of 
a = —1.85 ± 0.05 as our best determination. This is signifi- 
cantly (at the 3<r level) smaller than the "universal" initial 
CMF slopes of a — — 2 often found in interacting and star- 
burst galaxies, and used as the basis for theoretical models of 
the evolution of young cluster populations (see, e.g., de Grijs 
ct al. 2003d for a review). We will discuss the implications 
of this result in the context of our discussion of Fig. 8 below. 
We note that for log(M c i/M Q ) > 3.0, log(i dis /yr) > 9.3; in 
this age range, the clusters that would have been affected by 
ongoing (semi-instantaneous) disruption have already faded 
to below the adopted completeness limit, so that we con- 
clude that here we are indeed observing the initial CMF 
slope. We note that the low density of the LMC field, while 
to some small extent effective in tidal stripping and evapo- 
ration, does not lead to a significant breakdown of our as- 
sumption of (almost) instantaneous disruption: most of the 
evaporated and stripped stars will be of low mass, which 
hence contribute negligibly to the integrated luminosities 
we use for our analysis. 

Using the value of a = —1.85, and the fading parame- 
ter £ = 0.648 (where fading of a cluster's flux, Fx, is defined 
as Fa ~ t - '') from stellar population synthesis (Boutloukos 
& Lamers 2003), the predicted power-law slope for the low- 
mass range is (1/C) — \ct\- This slope is shown as the dash- 
dotted line in Fig. 6b, and matches the observed mass dis- 
tribution remarkably well (we have only applied a vertical 
shift to the slope in order to match the cluster numbers in 
the distributions). 

Alternatively, we can assess the CMF and its possi- 
ble evolution by examining the CMFs of well-defined and 
well-understood subsamples covering restricted mass ranges. 



Table 2. Mass range-dependent LMC cluster mass function 
slopes, based on magnitude-limited sampling (Fig. 6b). 



log(M c i/M ) mir 


CMF slope 


dAW 
d(M cl /M ) 


dlog(AT cl )/ 
dlog(M cl /M ) 


2.25 


-1.59 ± 0.08 


-0.59 ± 0.08 


2.50 


-1.67 ± 0.06 


-0.67 ± 0.06 


2.75 


-1.73 ± 0.06 


-0.73 ± 0.06 


3.00 


-1.88 ± 0.05 


-0.85 ± 0.05 


3.25 


-1.83 ± 0.05 


-0.83 ± 0.05 


3.50 


-1.85 ± 0.06 


-0.85 ± 0.06 


3.75 


-1.86 ± 0.08 


-0.86 ± 0.08 


4.00 


-1.84 ± 0.09 


-0.84 ± 0.09 


4.25 


-1.81 ± 0.12 


-0.81 ± 0.12 


4.50 


-1.80 ± 0.15 


-0.80 ± 0.15 


4.75 


-1.78 ± 0.21 


-0.78 ± 0.21 


5.00 


-1.55 ± 0.19 


-0.55 ± 0.19 



Following Fall et al. (2005), in Fig. 8 we present mass- 
limited LMC cluster subsamples, which are potentially phys- 
ically more informative than magnitude-limited subsamples. 
Mass-limited sample are less biased toward young clusters 
than magnitude-limited samples. This is a novel approach 
for the LMC star cluster system. In essence, a mass-limited 
statistically complete sample of clusters implies an imposed 
age limitation as well, so that for increasing masses, the up- 
per age limit increases, and so does the median age of the 
subsample. The mass and age ranges used for the construc- 
tion of the CMFs in Fig. 8 are listed in Table 3, as well as 
in the individual panels in the figure. 

The error bars indicate the simple Poissonian uncer- 
tainties, and the vertical dotted lines indicate the minimum 
mass limits adopted for a particular cluster subsample; in 
all cases the corresponding age limits are well below the ex- 
pected disruption time-scales for even the least massive clus- 
ters in a given sample. In Table 3 we also include the CMF 
slopes we derive from each cluster subsample. For samples 
up to log(Age/yr) < 8.7, we use a mass fitting range of 
log(M c i/M ) > 3.0; for the remaining subsamples we use 
the full mass range available for the fits. We note explicitly 
that using a lower age limit of log(Age/yr) = 6.8 or 7.2, in 
order to avoid the artifact caused by the limited age reso- 
lution of the isochrones used for the analysis does not alter 
these results significantly (or may even strengthen these re- 
sults, albeit marginally). 

The key result is that, for clusters with log(M c i/M Q ) > 
3.0 and log(Age/yr) > 7.80 (i.e., for all but the two youngest 
subsamples), the CMF slopes are very close to unity (see also 
H03 for similar results for the high-mass wings of their age- 
limited CMFs 1 ), which in this parameter space corresponds 
to an initial CMF slope of a ~ — 2 (see also the final column 
in Table 2 for a direct comparison) . However, the youngest 

1 We note that H03 did not take the age and mass-dependent 
completeness limits into account for the construction of their 
CMFs (their fig. 6). However, a close comparison between the 
mass fitting ranges they used on the one hand with Fig. 7 on the 
other reveals that they determined the CMF slopes well inside 
the statistically complete mass range for each age-limited cluster 
subsample. It is therefore not surprising, indeed rather encourag- 
ing, that the results we obtain for a larger number of subsamples 
for different age and mass ranges agree very well with those of 
H03. 
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Figure 8. CMFs for statistically complete LMC cluster subsamplcs. Age and mass ranges arc indicated in the panel legends; the vertical 
dotted lines indicate the lower mass limits adopted. Error bars represent simple Poissonian errors, while the best-fit lines were obtained 
for log(JW c i/MQ) > 3.0 for log(Age/yr) < 8.7, and for the full available mass range for older maximum ages. Table 3 summarises the 
best-fit parameters. 



two subsamples exhibit significantly shallower slopes for the 
same mass range (see also Elmegreen & Efremov 1997), 
while for log(M c i/M©) < 3.0 the observed CMFs turn down 
even more significantly below the linear slopes indicated. We 
note that the observed effect goes well beyond that expected 
from uncertainties in the data's completeness limit: an un- 
certainty in the completeness limit of 0.5 mag translates 
to an uncertainty in mass of only Alog(M c i/M Q ) = 0.2. 
As the Poissonian error bars indicate, this flattening of the 
CMF slope is unlikely to be caused by small-number statis- 
tics either, while cluster disruption or variable external tidal 
effects are expected to have had a more important effect on 
the low-mass extreme of older clusters than on the younger 
clusters for which we notice this discrepancy. Thus, this ob- 
servation is contrary to theoretical predictions if the initial 
CMF of the LMC cluster system were a power-law distribu- 
tion modified by cluster disruption processes. We recall that 



based on our magnitude-limited sample analysis (Fig. 6b), 
we also concluded that the initial CMF slope appeared to be 
significantly shallower, at the 3<r level, than the power-law 
slope, q = —2, expected for young star cluster systems. 

Since we observe this behaviour towards the low-mass 
end for all mass and age subsets of the full mass-limited 
cluster sample that include masses log(M c i/MQ) < 3, as 
well as in our magnitude-limited sample, random stochastic 
effects are unlikely to be the primary cause (see also Girardi 
& Bica 1993; Santos & Frogel 1997). It is likely that this 
is a real effect (see also Elmegreen & Efremov 1997), and 
that the CMF slopes flatten significantly for younger ages 
and lower-mass clusters in the LMC. Elmegreen & Efremov 
(1997) argue that the younger clusters are mostly unbound 
OB associations (supporting Bica et al. 1996), of which some 
90 per cent will disperse by the time their constituent stars 
will reach an age of ~ 10 8 yr. This is consistent with modern 
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Table 3. LMC cluster mass function slope as a function of min- 
imum mass and maximum age, based on mass-limited samples 
(Fig. 8). 



log(M cl /M Q ) min 


log(Age/yr) max 


[dlo 


CMP slope 
g(JV c i)/dlog(M c i/M )] 


2.2 


7.50 




-0.80 ± 0.10 


2.4 


7.80 




-0.68 ± 0.06 


2.6 


8.10 




-1.00 ± 0.10 


2.8 


8.40 




-0.94 ± 0.10 


3.0 


8.70 




-0.98 ± 0.08 


3.2 


8.95 




-0.95 ± 0.06 


3.4 


9.10 




-1.15 ± 0.11 


3.6 


9.35 




-1.18 ± 0.14 


3.8 


9.55 




-1.10 ± 0.12 
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ideas on the formation and dissolution of star clusters within 
the first ~ 10 7 yr of their existence; most (~ 70 — 90 per 
cent) of these newly formed clusters will disperse on these 
time-scales, a process coined "infant mortality" (e.g., Boily 
& Kroupa 2003; Vesperini & Zepf 2003; Whitmore 2004; 
Bastian et al. 2005; Mengel et al. 2005; see also Tremonti et 
al. 2001). If the process of infant mortality is mass dependent 
(but see Whitmore 2004 for counterarguments) , in the sense 
that the lowest-mass clusters will dissolve preferentially, the 
result will be a flattening of the CMF slopes with increasing 
mean age, which is contrary to the apparent change of slope 
in Fig. 8. Whether or not the youngest (< 10 Myr old) clus- 
ters will dissolve because of this infant mortality scenario, 
the significant differences among the LMC's CMF slopes as 
a function of age are predominantly driven by the youngest 
subset of our cluster sample. 

Thus, these results may imply that the initial CMF 
slope of the combined (i.e., both bound and unbound) LMC 
cluster system is not well represented by a power-law, al- 
though we cannot disentangle the unbound from the bound 
clusters at the youngest ages. In addition, we recently pre- 
sented observational evidence and theoretical arguments 
against an initial power-law CMF in M82's intermediate-age 
starburst region, M82 B (de Grijs, Parmentier & Lamers 
2005), and in favour of an initial log-normal CMF in the 
Antennae interacting system, NGC 4038/39 (Anders et al. 
2005). Our detailed analysis of the LMC cluster mass dis- 
tributions as a function of age and mass may therefore have 
uncovered supporting new evidence that star clusters - at 
least in the low-density environment of the LMC - may not 
form following a power-law distribution that holds strength 
down to masses much below a few xlO 3 M , where the 
power-law fits to the subsamples in Fig. 8 break down. 



5 SUMMARY AND CONCLUSIONS 

By combining integrated properties with resolved stellar 
population studies, the LMC cluster system offers the unique 
chance to independently check the accuracy of age (and 
corresponding mass) determinations based on broad-band 
SEDs. In this paper, we have reanalyzed the broad-band 
LMC cluster SEDs based on the data of Massey (2002) and 
H03, using a newly developed SED analysis approach. We 
compare our new age determinations with (i) those of H03 
using the same data set but a different approach, and (ii) 
those of Pietrzyriski & Udalski (2000) using CMD fitting, in 



order to set the tightest limits yet on the accuracy of (ab- 
solute) age determinations based on broad-band SEDs, and 
therefore on the usefulness of such an approach. 

We note a significant systematic effect between the age 
differences of H03 on the one hand, and those of both the 
OGLE-II team and our own redeterminations on the other. 
It appears that these systematic differences are caused by 
H03's conversions of the photometry to a different filter sys- 
tem. We emphasize and warn that the actual filter systems 
used for the observations should be used for the most ac- 
curate parameter analysis, instead of using filter conversion 
equations, in order to achieve more accurate derivations of 
the cluster ages and the corresponding masses. 

Based on this comparison, and additionally on a de- 
tailed assessment of the age-metallicity and age-extinction 
degeneracies, we conclude that our broad-band SED fits 
yield reliable ages, with statistical absolute uncertainties 
within Alog(Age/yr) ~ 0.4 overall. Thus, in addition to 
our conclusion in de Grijs et al. (2005) that we can re- 
trieve prominent features in the cluster age distribution to 
within A(log(Age/yr)) < 0.35 using a variety of approaches 
based on broad-band SEDs modelling, we have now also 
shown that the associated intrinsic statistical uncertainties 
involved in cluster age determinations are of a very similar 
magnitude. 

The LMC's CFR has been roughly constant outside of 
the well-known age gap between ~ 3 and 13 Gyr, when the 
CFR was a factor of ~ 5 lower (assuming a roughly con- 
stant rate during the entire period) . There are no clear ob- 
servational signatures of an enhanced CFR associated with 
the last tidal encounter between the LMC and the SMC, 
while we argue that the combination of the relevant time- 
scales, i.e., the LMC's rotation period and the time since 
the last LMC-SMC encouter, has been insufficient to wash 
out any such signatures, if they had been present. An alter- 
native triggering mechanism for the young(er) clusters may 
be needed. 

Using a simple approach to derive the characteristic 
cluster disruption time-scale, we find that log(£4 ls /yr) = 
9.9 ± 0.1, where t dis = tf s (M cl /lO 4 M ) ' 62 , for the LMC 
cluster system. This is consistent with earlier, preliminary 
work for a smaller cluster sample. This long characteristic 
disruption time-scale implies that hardly any of our LMC 
sample clusters are affected by significant disruptive pro- 
cesses, so that we are in fact observing the initial CMF. 
Using a variety of complementary techniques, we conclude 
that the older cluster (sub) samples show CMF slopes that 
are fully consistent with the a ~ — 2 slopes generally ob- 
served in young star cluster systems. The youngest clusters 
in our sample show shallower slopes, at least below masses 
of a few xlO 3 Mq, and possibly evidence for a turn-over. 
This is contrary to dynamical expectations and may imply 
that the initial CMF slope of the LMC cluster system as a 
whole is not well represented by a power-law down to the 
lowest masses, although we cannot disentangle the unbound 
from the bound clusters at the youngest ages. 
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